Universal decoherence induced by an environmental quantum phase transition 
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Decoherence induced by coupling a system with an environment may display universal features. 
Here we demostrate that when the coupling to the system drives a quantum phase transition in the 
environment, the temporal decay of quantum coherences in the system is Gaussian with a width 
independent of the system-environment coupling strength. The existence of this effect opens the way 
for a new type of quantum simulation algorithm, where a single qubit is used to detect a quantum 
phase transition. We discuss possible implementations of such algorithm and we relate our results 
to available data on universal decoherence in NMR echo experiments. 

PACS numbers: 



The coupling between a quantum system and its en- 
vironment leads to decoherence, the process by which 
quantum information is degraded. Decoherence plays a 
crucial role in the understanding of the quantum to classi- 
cal transition |l|. It also has practical importance: its un- 
derstanding is essential in technologies that actively use 
quantum coherence, such as quantum information pro- 
cessing 0]. In general, the timescale tdec of decoherence 
depends on the system-environment coupling strength, 
which we arbitrarily denote A. For example, in the well 
studied case of quantum Brownian motion (where the en- 
vironment consists of a large number of non-interacting 
harmonic oscillators), quantum coherence generally de- 
cays exponentially with a rate 1/idec proportional to A^ 
[2j. In this letter we describe a class of systems with a 
drastically different behavior: Gaussian decay of coher- 
ence with a rate independent of A. This independence 
signals a universal behavior whose study is the aim of 
this work. In general, one should avoid building physical 
quantum information processing devices in presence of 
universal decoherence. However, we show that universal- 
ity is a powerful property we can use to our advantadge: 
by detecting decoherence in the universal regime we can 
extract valuable information about the environment. 

Environment-independent decoherence rates are also 
found in other circumstances. For example, systems 
with a classically chaotic Hamiltonian display a "Lya- 
punov regime" where the decay is exponential and given 
by the Lyapunov exponent of the underlying classical 
dynamics |j, |5j- These models are also often used to 
represent a complex environment. In fact, chaoticity is 
the widespread explanation |3, la| for the perturbation- 
independent decay of polarization detected in recent 
NMR echo experiments J3 (where, however, a non- 
exponential but Gaussian decay is actually observed). 
Our findings are different from the usual exponential 
Lyapunov regime: we discuss systems where the univer- 
sal (independent of A) decoherence is Gaussian. In our 
model, the complexity and sensitivity of the environment 



arise from the susceptibility of the environmental spec- 
trum to the system's state. The relation between our 
results and the experiments of Ref. [3| will also be dis- 
cussed below. 

Let us consider a spin 1/2 particle (a qubit) coupled 
to an environment that is "structurally unstable" with 
respect to the system state (in a sense that will be made 
clear below). The model we discuss is a generalization 
of the one studied by Quan et al [8j, who showed that 
an environment at the critical point of a quantum phase 
transition is highly efficient in producing decoherence. 
Below, we will not only generalize the results of [Sj but 
also show that in these circumstances universal decoher- 
ence arises naturally. We assume that the system and 
the environment evolve under the Hamiltonian 

Use - Is ® Hb + |0) (0| ® Hao + |1) (1| ® Hai- (1) 

Here, the operators TLe: T^Xo ^-iid Tt\^ act on the Hilbert 
space of the environment. If the system is in state 
\j) {J — 0; l)j the environment evolves with an effec- 
tive Hamiltonian Tij — He + "Hx {Xj is the system- 
environment coupling strength). Considering the initial 
state |*5£(0)) = (a |0)+6 |1)) |£:(0)) , the evolved reduced 
density matrix of the system is 



Ps{t) 



= Tr£|*5£W)(*5£(i)l 

= |anO)(0|+a6Mt)|0)(l| 

+ a*6r*(t) |l)(0| + |5ni)(l| 



(2) 



The off-diagonal terms of this operator are modulated 
by the decoherence factor r(t): the overlap between two 
states of the environment obtained by evolving the ini- 
tial state |i?(0)) with two different Hamiltonians, i.e. 
r{t) = (£(0)1 e*^o*e~''"i* 15(0)). Moreover, assuming 
that the initial state of the environment is the ground 
state l^o) of Wo [^^ the decoherence factor r{t) is, up 
to an irrelevant phase factor, identical to the so-called 
survival probability amplitude 



r{t) = (50 1 e" 



'iHit 



.9o). 



(3) 



Let us first analyze models where both Hamiltonians 
Tij (j = 0, 1) can be diagonaUzed in terms of a suitable 
set of fermionic creation and annihilation operators 7^ : 



n, 



N 
k=l 



l^'^l^' 



(4) 



Furthermore, we assume that the operators appearing in 
the two Hamiltonians Jij can be connected by a Bogoli- 
ubov transformation of the form 



1^ 



I \ W ■ ■ I \ (o)t 
cos(afc)7;: -is\n{ak)X-k ' 



(5) 



where the angles at define the Bogoliugov coefficients. 
Notice that this expression only includes mixing between 
modes with opposite values of the index k. Our treatment 
can be extended to more complicated situations, but we 
limit first to the simplest non-trivial case, where it is 
possible to relate the ground states \g) , of Tij as 

\g), = n [^cos(a,) + sinK)7I')V'i^ 1.9) 1 • (6) 

/c>0 



Under these assumptions the decoherence factor is 



r{t) = n (cos'(afc) 



sin^(afe)e" 



(7) 



fc>0 



Surprisingly, r(i) is completely analogous to the one 
found when studying non-interacting spin environments 
[lOj. In that case, the index k labels the different envi- 
ronmental spins and the corresponding Bogoliubov coef- 
ficients define their initial states. 

Under reasonable assumptions on the angles ak and 
the energies e,, , we can go further and - using the ideas 
developed in [iflj - obtain a simple form for the temporal 
evolution of the overlap r{t). To illustrate our procedure, 
let us analyze first an oversimplified case: suppose that 
the energies of all the modes are the same, i.e. t]. — t. 
In the simplest case ak = 7r/4, the overlap oscillates as 
r{t) = (coset)^/^. The same result is recovered as a con- 
secuence of the law of large numbers if the angles au are 
spread over the entire circle. In fact, \r{t)\'^ ~ jcosei 
if the following Lindenberg conditions are satisfied 
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i- ^ cos^ ak ^ 1/2 

k 

= ^sin^ 2a. (.«)%> 



(8) 



The first condition is satisfied when the angles are ran- 
domly distributed. The second one imposes a finite vari- 
ance for the "quantum walk" in which a step of length 
+f-k (— Efe) is taken with probability cos^ a. (sin a.). 
e, the condition takes the form s?, ^ 1, 



When 



.(1) 



and it is met when there is a sufficiently large number of 
modes for which sin 2ak does not vanish. 

A more realistic situation is when the energies ej^, take 
values in a given spectral band. When the energies are 
distributed with a vanishing mean value, the decay of 
r(i) is Gaussian with a width given by s^ defined in 
© [ifl • Consider the more general case where the ener- 
gies are distributed about an arbitrary mean value, i.e. 
Ej, = e + 5k (where 5k has zero mean). We now define 
the dispersion s]^ as the cumulative variance of the fluc- 
tuations of the energy, i.e. 5^ = X^fc ^i^^ 2afc J^. We 
find that, in general, when conditions (jS)) hold (replac- 
ing s^ by s\), r(t) is described by a Gaussian envelope 
modulating an oscillating term. 



\r{t)\'' ^ eM-sli t')\cos{et)f'^ . 



(9) 



In general, when the operators 7^ and 7^, are simi- 
lar, the angles a. are small and © do not hold: there is 
almost no decoherence. However, a drastic difference in 
the nature of the eigenstates of TLq and Tii can only be ac- 
counted for with ak varying in the full range [0, 27r). This 
occurs when the environment suffers a quantum phase 
transition when A is varied. Thus, denoting Ac the criti- 
cal point of the transition, for Aq ^ Ac -C Ai we expect 
the decoherence factor to behave as indicated in 0. In 
many cases, s^ is only given by the properties of the 
environment Hamiltonian, and thus the decay of r(i) be- 
comes universal (independent of A). 

An important model encompassed by assumptions Q) 
and (jSl is an Ising chain transversely coupled to a central 
spin 8] (which plays the role of the system). In this case 
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-•/(E'^^^^i-^^E^?)- 



(10) 



The Bogoliubov coefficients and the energies are [llj. 



2 J Jl + \] - 2Xj cos(27rfc/7V) (11) 



^^' 



2ak = (0fe(Ai)-0fe(Ao)), 



(12) 



N 



where the angles ^.(A) are defined from tan(6'fe) — 
sm{2Trk/N)/{X - cos{2Trk/N)). In this model Ac = 1. 

When Ai 3> 1 and Ao < 1, the angles afc(A) w 
wk/N and the Bogoliubov coefficients satisfy conditions 
0. Moreover, the energies e^ ' are distributed between 
|Ai — 1| and |Ai -I- 1|, which gives s^ ~ N. Therefore, the 
width of the Gaussian envelope is independent of Xi. 

In Fig. 1 we display r{t) for the case Aq = 0, showing 
the accuracy of Eq. ©. The universality of the envelope 
is a clear indication of the quantum phase transition. 
However, the oscillations (whose frequency depends on 
Ai) are not universal. Yet, it is possible to eliminate 
them by performing a spin-echo experiment: first, evolve 
the system coupled to the Ising chain environment for a 
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FIG. 1: Decoherence factor for the Ising model Hamiltonian 
as the environment with 50 spins for (a) Ai = 2, (b) Ai = 5, 
(c) Ai = 10, and (d) Ai = 40. In (a) the two top curves 
included for comparison have Ai = 0.1 and Ai = 0.5. In all 
plots Ao = 0, and the dashed line is the predicted universal 
Gaussian envelope. On the right, the same values of Ai for 
recho{t), which eliminates the Ai dependent oscillations. 



time t. At this time, flip the environmental spins in the x 
direction (e.g. with an rf-pulse that applies a 7r-rotation 
around the z-axis). Finally, evolve for another time t. 
The total evolution of the environment can be described 
by using the Hamiltonian Tii = He + Hai from time 
to t, and Ti-i = He — "Hx^ from time t to 2t. Thus, in 
this echo experiment the decoherence factor is given by 



rechoC^t) = {gl^e' 



-iH- 



-iHit 



9)0 ■ 



(13) 



This overlap is simply computed using the Bogoliubov 
transformation that connect the modes diagonalizing 
the Hamiltonians 7i_i and Ho- If we denote 7^. 
the modes of H-i, the Bogoliubov coefficients associ- 
ated with the corresponding angles ak are such that 
7^ = cos((5fe)7^, — zsin(afe)7_^. . The analytic form 
for the overlap Vechoit) is simplified introducing the sum 



Bogoliubov angles, a^, — fifc ± ctfc. We obtain 



^echoy-^^) 
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k>0 



cose^ t cos^q:[, +cose^ 't sin'^ a^ 



4^'^ 
4-'^ 



4-^' 
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4^^ 



(14) 



For the case of the Ising model the expression can be eval- 
uated explicitely. In the limit of large values of A, one 
can obtain an approximate behavior using similar argu- 
ments as above Il2lj . Thus, the dominant contribution to 
the echo-overlap is 



recho{t) « exp(-sjvi ) 1 - 



m 

A 



sin(Ai) 



(15) 



where ii'(i) = 2J2kSHA~^t)cos{2Trk/N)sm'^{2TTk/N). 
In Fig. 1 we show how the accuracy of this expression 
increases with A. 

To test the generality of our results against the restric- 
tiveness and uncontroUability of assumptions Q and Q , 
we study a system in the opposite end of the spectrum: 
the Bose- Hubbard model (BHM) [11|, with Hamiltonian 



Hbh = -g 



ala-i 






alttnialan - 1). (16) 



<»j> 



and difference of the energies, e 



(±) 



e^'^±e^"'\andthe 



Here a„ are boson anihilation operators in site n of a 
discrete lattice. For g :^ u, the system behaves as a 
superfluid of non-interacting particles. In the opposite 
regime, u ^ g, the interaction term dominates and the 
ground state is Mott-insulator like. This model cannot 
be cast in terms of fermionic operators as in Q, in fact, 
no analytic solution is known. Furthermore, the bosonic 
nature of the particles also conflicts with (O- The BHM 
has practical relevance because it can be experimentally 
simulated using cold neutral atoms in an optical lattice 
|l3| . We calculate r{t) numerically for a spin 1/2 coupled 
to the hopping term of the BH Hamiltonian, that is, we 
take g = X. In Fig. 2 we show the decoherence factor 
for several values of A for a BHM with a fixed number of 
bosons. The same overall behavior of the Ising chain is 
observed: a universal Gaussian envelope (independent of 
A) modulating an oscillation with frequency proportional 
to A. The very different nature of the BHM hints at a 
more general validity of our results. 

A Gaussian decay of coherence with a rate independent 
of the coupling to the environment was indeed observed 
in NMR polarization echo experiments [3- Arguing 
on the complexity of the experimental many-body sys- 
tem, these results have been related to the environment- 
independent decoherence predicted in classically chaotic 
Hamiltonians |j, |^ |y|. The experimental situation is 
quite different from the one we considered here: the de- 
coherence factor is measured after an echo created by a 




t [h/u] 



FIG. 2; Decoherence factor for the Bose-Hubbard model 
Hamiltonian as the environment, with 6 particles in a lat- 
tice of 6 sites, for (a) J = 5, (b) J = 10, (c) J = 20, and (d) 
J = 50. In all plots the dashed line is the universal Gaussian 
envelope whose width is numerically obtained. 



change of sign of the environment Hamiltonian, and not 
the system-bath interaction. Our model points to a dif- 
ferent way of introducing complexity and sensitivity in 
the environment: a quantum phase transition. Further 
research using this approach might explore more realistic 
models that account for all the details of the experiments. 

The universal decoherence regime of this work can 
also be understood using analogies to the regime of 
strong perturbations of the survival probability, Eq. ^. 
Indeed, r(i) is the Fourier transform of the strength 
function or local density of states (LDOS), L{E) — 
En I (Sol^n) \'^5{E - En), where |(/)„) are the eigenvec- 
tors of Til and £'„ its eigenenergies. In typical LDOS 
studies. Til differs from Tio by a perturbation. In com- 
plex systems (e.g. random matrices, or classically chaotic 
Hamiltonians) for sufficiently strong perturbations jgo) 
is a random superposition of the |0„) states. Therefore, 
the LDOS becomes independent of the perturbation: it 
equals the full density of states of Tti . In our model, the 
saturation of the LDOS when TYq and Tii are on both 
sides of the quantum phase transition occurs because of 
the radically different nature of the eigenstates. In con- 
trast to our results, Refs. [lj| have found that complex 
systems give an LDOS with a Lorentzian shape, leading 
to an exponential decay of r{t). 

Universal decoherence can be harmful for quantum in- 
formation applications. However, it can be a useful tool 
to extract information about a critical system, e.g. its 
spectral structure or the critical point of its quantum 
phase transition. The latter example can be thought 
of as a "critical point finding" algorithm in a one-qubit 



quantum computer: in systems where the spectrum is 
not shifted by the coupling (which gives the oscillatory 
cos{Xt)^ term), the critical point can be simply obtained 
as the A value for which one observes the onset of uni- 
versality. Otherwise, the oscillation term obscures the 
critical point. In these cases one can instead couple the 
system weakly to the environment, and drive the tran- 
sition with an external parameter (as in Ref. [Sj). The 
critical point is then signaled by the A value for which 
there is a maximum decoherence decay. A demonstration 
of this algorithm can be performed in an NMR setting 
simulating the Ising Hamiltonian studied above [l5j. 

We have shown that when the coupling to the sys- 
tem drives a quantum phase transition in the environ- 
ment, the decoherence factor decays as a Gaussian with 
an environment-independent width. We showed numer- 
ically that our findings are more general than what can 
be expected from the analytical approximations we used. 
Our results could lead to an alternative interpretation of 
hitherto unexplained NMR experimental results on en- 
vironment independent decoherence rates. Finally, we 
discussed how the universal behavior of the decoherence 
factor can be used to study critical systems in a novel 
simulation algorithm for one-qubit quantum computers. 
We acknowledge fruitful discussions with W.H. Zurek. 
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